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ABSTRACT 

Retinoic acid-related orphan receptors (RORs) and 
the basic helix-loop-helix-PAS transcription factor 
Npas2 have been implicated in the control of circa- 
dian rhythm. In this study, we demonstrate that 
RORy directly regulates Npas2 expression in vivo. 
Although the rhythmicity of Npas2 mRNA expression 
was maintained in RORy~'~ mice, the peak level 
of expression was significantly reduced in several 
tissues, while loss of RORa had little effect. 
Inversely, overexpression of RORy in hepatoma 
Hepa1-6 cells greatly induced the expression of 
Npas2. RORy-activated Npas2 transcription directly 
by binding two ROREs in its proximal promoter. 
ChIP analysis demonstrated that RORy was 
recruited to this promoter in the liver of wild-type 
mice, but not RORy-deficient mice. Activation of 
Npas2 correlated positively with chromatin acces- 
sibility and level of H3K9 acetylation. The activation 
of Npas2 by RORy was repressed by co-expression 
with Rev-Erba or addition of the ROR inverse 
agonist T0901317. Npas2 expression was also 
significantly enhanced during brown adipose differ- 
entiation and that this induction was greatly sup- 
pressed in adipose cells lacking RORy. Our results 
indicate that RORy and Rev-Erba are part of a 
feed-back loop that regulates the circadian expres- 
sion of Npas2 suggesting a regulatory role for 
these receptors in Npas2-dependent physiological 
processes. 

INTRODUCTION 

The retinoic acid-related orphan receptor (ROR) subgroup 
of nuclear receptors consists of three members, RORa, -P 
and -y (NR1F1-3 or RORA-C) (1). As a result of 



alternative promoter usage and/or exon splicing each 
ROR gene generates several isoforms that differ only in 
their amino terminus. These variants exhibit a distinct 
pattern of tissue-specific expression and are involved in 
the regulation of different physiological processes and 
target genes. The RORs bind as a monomer to ROR 
response elements (ROREs) in the promoter regulatory 
region of target genes (2,3). Transcriptional regulation by 
RORs is mediated by the recruitment of co-repressor and 
co-activator complexes which, through their his- 
tone(de)acetylase activities, induce changes in chromatin 
conformation (1,4). Recent studies have demonstrated 
that RORs function as ligand-dependent transcription 
factors. A\l-trans retinoic acid and several other retinoids 
act as partial antagonists of RORP and T0901317 and 
7a-hydroxycholesterol as inverse agonists of RORa and 
RORy, while cholesterol sulfate has been reported to 
function as an RORa agonist (5-8). 

RORs regulate several important physiological 
processes and have been implicated in a number of 
pathologies (1). RORa is critical for cerebellar develop- 
ment and bone formation (9,10), while RORP regulates 
functions in the brain and retina (11,12). RORy plays a 
key role in lymph node development and thymopoiesis 
(13,14). Furthermore, both RORa and RORy are 
involved in regulating various metabolic pathways, in- 
flammatory responses and immune functions, including 
Thl7 cell differentiation (1,15-19). 

In addition to these functions, all three RORs have been 
implicated in the regulation of various aspects of circadian 
rhythm (1,11,20-27). Circadian rhythm is fundamental in 
the regulation of a wide variety of physiological and 
behavioral activities. In mammals, the central clock in 
the suprachiasmatic nucleus (SCN) integrates light-dark 
cycle input and synchronizes the autonomous oscillators 
in peripheral tissues (28-30). The molecular basis of the 
circadian clock involves interlocking transcriptional/ 
translational feedback loops that regulate the rhythmic 
expression and activity of a set of core clock genes, 
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including the basic helix-loop-helix/Per-Arnt-Sim (PAS) 
domain-type transcriptional activators brain and muscle 
aryl hydrocarbon receptor nuclear translocator 
(ARNT)-like {Bmall), circadian locomotor output cycles 
kaput (Clock), two cryptochrome (Cry) and three period 
(Per) proteins. The positive loop consists of the basic 
helix-loop-helix/PAS-type transcription activators 
Bmall and Clock, while Cry and Per are involved in the 
negative control of the oscillator (31-36). The repression 
of Bmall is mediated by the nuclear orphan receptors 
Rev-ErbAa and Rev-Erbp, which function as transcrip- 
tional repressors, while RORs have been implicated in 
the positive regulation of several clock genes, including 
Bmall and Rev-Erba (1,22-24,26,37,38). The circadian 
clock subsequently regulates the rhythmic expression of 
downstream genes involved in timing the rhythmic regu- 
lation of numerous functions, including behavior, repro- 
ductive and neuroendocrine functions, and metabolism. 

In this study, we identify neuronal PAS domain protein 
2 (Npas2), a paralogs of Clock, as a new ROR target gene. 
Npas2 heterodimerizes with Bmall and regulates the tran- 
scription of circadian genes (39-41). Npas2 and Clock 
have been reported to exhibit overlapping functions. 
Single nucleotide polymorphisms (SNPs) in Npas2 have 
been linked to increased risk of cancer, unipolar major 
depression, metabolic syndrome and hypertension (42- 
47). In this study, we demonstrate that the loss of particu- 
larly RORy significantly reduced the peak level of Npas2 
mRNA expression in several tissues without affecting the 
rhythmicity of Npas2 expression. We further show that 
RORy regulates Npas2 transcription directly through 
two ROREs in its proximal promoter. This activation is 
inhibited by the inverse RORy agonist T0901317 and 
by Rev-Erba. Collectively, these data indicate that both 
RORy and Rev-Erba are part of the mechanism that regu- 
lates the oscillatory expression of Npas2 and suggest a 
regulatory role for these receptors in Npas2-dependent 
physiological processes and pathologies, including circa- 
dian behavior disorders, metabolic syndrome and 
tumorigenesis. 



MATERIALS AND METHODS 

Experimental animals 

Heterozygous C57BL/6 staggerer (RORct +/sg ) mice were 
purchased from Jackson Laboratories (Bar Harbor, ME, 
USA). The generation and initial characterization of 
C57BL/6 RORy~ l ~ mice were described earlier (13). 
Mice were supplied ad libitum with NIH-A31 formula 
and water as described earlier (16) and maintained at 
25° C on a constant 12 h light: 12 h dark cycle. Littermate 
wild-type (WT) mice were used as controls for both 
ROR-deficient models. For cold-induced thermogenesis 
mice were placed at 4°C for 4h. All animal protocols 
followed the guidelines outlined by the NIH Guide for 
the Care and Use of Laboratory Animals and were 
approved by the Institutional Animal Care and Use 
Committee at the NIEHS. 



RNA isolation 

To study gene expression during circadian time, tissues 
were excised from WT, RORy~ ! ~, and RORa sg/sg mice 
every 4h over a period of 24 h, processed overnight in 
RN Alater® solution (Ambion, Austin, TX, USA) at 
4°C, and then stored at — 80°C until use. Tissues were 
then homogenized in RLT buffer (Qiagen, Valencia, CA, 
USA) with a Polytron PT-3000 (Brinkmann Instruments, 
Westbury, NY, USA). RNA was then extracted using a 
QIAshredder column and RNeasy Mini kit (Qiagen) ac- 
cording to the manufacturer's instructions. Hepal-6 cells 
and brown adipocytes were lysed in lysis (RLT) buffer and 
RNA extracted as described for tissues. The brown adipo- 
cytes were treated for 24 h (from differentiation Day 4 to 
5) with and without 10 uM of the inverse ROR agonist 
T0901317 (7) (Sigma-Aldrich, St Louis, MO, USA) before 
RNA extraction. 

Quantitative reverse transcription-polymerase chain 
reaction (QRT-PCR) analysis 

QRT-PCR analysis was performed using SYBR Green I 
or TaqMan system (Applied Biosystems, Foster City, CA, 
USA) to quantify gene expression. All primer and probe 
sequences are listed in Supplementary Table SI. The RNA 
was reverse-transcribed using a High-Capacity cDNA 
Archive Kit (Applied Biosystems). The reactions were 
carried out in triplicate using 20 ng of cDNA and the fol- 
lowing conditions: 2min at 45°C and lOmin at 95°C, 
followed by 40 cycles of 15s at 95°C and 60s at 60°C. 
All the results were normalized by the amount of 18S 
rRNA or Gapdh mRNA. 

Isolation of brown preadipocytes 

Brown fat preadipocytes were isolated from postnatal 
day 2 (PND2) WT and RORy~ l ~ brown adipose tissue 
by collagenase digestion as described earlier (48). 
Preadipocytes were immortalized by infection with pBabe 
retrovirus encoding SV40-Large T antigen. After selection 
with puromycin (2mg/ml), resistant clones were isolated. 
Several clones of the each genotype, WT, RORy +l ~ and 
RORy~ ~ were obtained. Cells were grown in DMEM 
supplemented with 10% fetal bovine serum (FBS), 1 nM 
thyroid hormone (T3) and 20 nM insulin (differentiation 
medium). Differentiation in confluent cultures was 
induced by the addition of 1 uM dexamethasone, 
0.125 mM indomethacin and 0.5 mM isobutylmethyl- 
xanthine. After 2 days incubation, the medium was 
changed to the normal differentiation medium and 
replaced every other day as described earlier (49). 

Hepal-6 and brown preadipocyte stable cell lines 

Hepal-6 cell lines stably expressing Flag-RORa4, 
Flag-RORyl and Flag-Rev-Erba were generated after 
transfection with the respective pLXIN plasmid DNA 
and subsequent selection in G418 (Clontech, Palo Alto, 
CA, USA). In the remaining of the article, RORa4 and 
RORyl are commonly referred to as RORa and RORy. 
Brown adipocyte cell lines BAT(E), BAT(RORy) and 
BAT(RORyE502Q) stably expressing empty vector. 
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RORy or the RORyE502Q mutant containing an inactive 
activation domain, respectively, were generated in the 
same way from immortalized preadipocytes. The expres- 
sion and nuclear localization of RORy protein were con- 
firmed by western blot analysis and immunochemical 
staining, respectively. Cell lines established from five indi- 
vidual clones were examined in QRT-PCR and chromatin 
immunoprecipitation (ChIP) analysis. 

ChIP assay 

ChIP assay was performed using a ChIP assay kit from 
Millipore (Billerica, MA, USA) according to the manufac- 
turer's protocol with a few minor modifications. Briefly, 
the livers isolated from WT, RORa sg/sg and RORy~'~ 
mice were homogenized with a polytron PT 3000 
(Brinkmann Instruments) and cross-linked by 1 % formal- 
dehyde for 20min at room temperature. After a wash in 
PBS, an aliquot of the cross-linked chromatin was 
sonicated and incubated overnight with anti-RORa or 
anti-RORy antibody (sc-6062 sc-28559, respectively; 
Santa Cruz Biotechnology, Santa Cruz, CA, USA). To 
analyze histone H3 lysine 9 (H3K9) acetylation within 
the Npas2 promoter region, ChIP assays were performed 
with an anti-H3K9Ac antibody (ab4441; Abeam, 
Cambridge, MA, USA). After incubation with protein G 
agarose beads for 2h, and several washes, DNA-protein 
complexes were eluted. The cross-links were reversed by 
overnight incubation at 65°C in the presence of 25 mM 
NaCl and subsequent digestion with RNase A and pro- 
teinase K, followed by purification of the chiped-DNA. 
The fold amount of the chiped-DNA relative to each 
input DNA was determined by QPCR. All QPCR reac- 
tions were carried out in triplicate. The procedure for 
ChIP analysis using Hepal-6 cells or brown adipocytes 
was similar as described for tissues except that 2 x 10 6 
cells were cross-linked with 4% formaldehyde for lOmin 
and the immunoprecipitation was performed with 
anti-Flag M2 affinity gel (Sigma-Aldrich). The sequences 
of the primers used for ChlP-QPCR are listed in 
Supplementary Table SI. 

Formaldehyde-assisted isolation of regulatory elements 

Formaldehyde-assisted isolation of regulatory elements 
(FAIRE) was performed as described earlier (50). 
Briefly, mouse livers were cross-linked and sonicated as 
described for the ChIP assay. Samples were centrifuged 
and the DNA in the supernatants isolated by three con- 
secutive extractions with phenol-chloroform-isoamyl 
alcohol (25:24:1). After the final extraction, the FAIRE 
samples were reverse-cross-linked as described for the 
ChIP assay. The enrichment of fragmented genomic 
DNA in the FAIRE samples relative to each input 
DNA was measured in triplicate by QPCR. 

Reporter gene assay 

To generate pGL4.10-A^av2(-1534/+81), the region of 
— 1534/+81 bp of the murine Npas2 promoter was cloned 
into the promoter-less reporter plasmid pGL4.10 
(Promega, Madison, WI, USA). The mutant pGL4.10- 
Ayms2(-1534/+81) plasmids, in which one of the two 



ROREs or both ROREs were mutated from AGGTCA 
to AAATCA, were generated using a Quickchange site- 
directed mutagenesis kit (Stratagene, La Jolla, CA, USA). 
The sequences were verified by DNA sequencing. Huh-7 
and HEK293T cells were co-transfected with pGL4.10- 
Npas2(-l534/+8\l pGL4.10, or pCMVP-Gal reporter 
plasmids and p3xFlag-CMV10-RORy, p3xFlag-CMV10- 
RORa or p3xFlag-CMV10-Rev-Erba expression plasmids 
as indicated using Lipofectamine 2000 (Invitrogen, 
Carlsbad, CA, USA). After 24 h incubation, the luciferase 
and P-galactosidase activities were measured by Luciferase 
Assay Substrate (Promega) and Luminescent P-galacto- 
sidase Detection Kit II (Clontech). All transfections were 
performed in triplicate, and each experiment was repeated 
at least twice. 



RESULTS 

Regulation of circadian expression of Npas2 by RORs 

Several studies have provided evidence for a role of RORs 
in the regulation of several circadian clock genes 
(1,11,20,23,24,26,27,38). In this study, we focused on the 
regulation of Npas2 by RORa and RORy. First, we 
analyzed the effect of the loss of RORoc or RORy expres- 
sion on the circadian expression of Npas2 in the liver, 
brown adipose tissue (BAT), kidney and small intestines 
using RORa g/sg and RORy~ h "mice. Consistent with a 
previous report, Npas2 mRNA showed a strong oscilla- 
tory pattern of expression in the liver of WT mice with 
peak expression between CT20 and CT0 and a low point 
at CT12 (Figure 1A and B) (51). However, the oscillatory 
expression of Npas2 was significantly altered in RORy~ l ~ 
liver. At CT0, when Npas2 is highly expressed in WT liver, 
expression of Npas2 mRNA was reduced by about 60% in 
the liver of RORy~ ! ~ mice (Figure 1A). Npas2 expression 
exhibited a similar robust oscillatory pattern in BAT, 
kidney, small intestines and white adipose tissue (WAT) 
as in WT liver (Figure 1C-H). However, between CT20 
and CT4 Npas2 expression was greatly decreased in BAT 
of RORy~~ mice (Figure 1C), reduced modestly in 
RORy~ l ~ kidney (Figure IE) and was slightly enhanced 
in the small intestines (Figure 1G). Little change was 
observed in the pattern of Npas2 expression in the liver 
and WAT from RORa g/sg mice (Figure IB and H), while 
a small shift in peak Npas2 expression was seen in 
RORa sg/sg kidney (Figure IF) and an increase in 
RORa xg/sg BAT. RORa is not expressed in small 
intestines. 

Earlier studies identified a certain degree of redundancy 
between the regulation of gene expression by RORa and 
RORy (16). To examine whether there was any func- 
tional redundancy between RORa and RORy, Npas2 
expression was examined in RORa sglsg RORy~ l ~ double 
knockout (DKO) mice. As shown in Figure II, at CT20 
Npas2 expression was reduced to a very similar degree in 
the liver, BAT and kidney of DKO mice as observed in 
RORy~ l ~ mice. Together, these results indicate that the 
circadian expression of Npas2 is largely regulated by 
RORy and to a minor extent by RORa and that this regu- 
lation exhibits a certain tissue specificity. 
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Figure 1. Comparison of the circadian expression of Npas2 in several peripheral tissues from WT, RORy'^ and RORa sglsg mice. (A-H) Various 
tissues from WT, RORy^ 1 ^ and RORct sg/sg mice (n = 4) were isolated every 4h over a period of 24 h and expression of Npas2 was analyzed by 
QRT-PCR. The 24-h expression pattern was double plotted. (I) Tissues were isolated from RORy^ 1 ^, RORa sg,sg RORy^ h DKO mice (n — 4) at 
CT20. Npas2 expression in DKO and RORy - ' - mice was normalized to the expression in littermate WT mice. (J) Circadian pattern of RORyl, 
RORal and RORa4 mRNA expression in liver from WT, RORy~^~ and RORa sg ' sg mice and Rev-Erba expression in WT mice (n = 4). Data present 
mean ± SD, */><0.05, **P<0.01, ***/>< 0.001 by ANOVA. 
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To be able to relate the oscillatory pattern of expression 
of Npas2 in liver to those of RORa and RORy, we 
examined their circadian pattern of expression. Figure 1J 
shows that particularly RORyl, the only RORy isoform 
expressed in liver, showed a robust oscillatory pattern of 
expression, while RORal, which is expressed at low levels, 
has a moderate and RORa4, which in liver is the predom- 
inant a isoform, exhibited a weak oscillatory pattern con- 
sistent with earlier observations (16,20,52). Loss of RORy 
or RORa expression had little effect on the circadian ex- 
pression of RORal or c/4 and RORyl mRNA, respective- 
ly, suggesting that RORa and RORy are regulated 
independently from each other (Figure 1J). The phase of 
the oscillatory expression of RORyl in liver (zenith at 
CT18) was slightly earlier than that of Npas2 (zenith at 
CT20) and the opposite of that of RORal and Rev-Erba 
expression (zenith at CT8) (Figure 1J). 

Circadian expression of Clock gene in RORa. sg/sg and 
RORy~l~ mice 

Npas2 and Clock are paralogues that share several func- 
tions (39 — 41). Each protein can form a heterodimer with 
Bmall. Moreover, Npas2 is able to functionally substitute 
for Clock in the master clock in the SCN and overlapping 
functions have been reported in the liver. A putative 
RORE has been identified in the human and mouse 
promoter region of the Clock gene, while EMSA 
analysis showed that RORs were able to bind to this 
RORE (27). To further examine the potential role of 
RORs in the regulation of Clock, we next determined 
the effect of the loss of RORa or RORy on its expression. 
In the liver, BAT, kidney and intestines Clock exhibited a 
less-robust oscillatory pattern of expression than Npas2 
with maximum expression around CTO and a similar 
circadian phase as Npas2 (Supplementary Figure SI). 
No significant differences in the circadian expression 
pattern of Clock were observed in the liver, BAT and 
WAT between WT and RORa sg/sg or RORy~ ! ~ mice 
(Supplementary Figure S1A-D and H). Interestingly, at 
CT20 and CTO Clock expression was greatly down- 
regulated in the kidney of RORa sglsg mice compared 
with WT kidneys (Supplementary Figure S1F), while 
loss of RORy significantly reduced the amplitude of the 
oscillatory expression of Clock in kidney (Supplementary 
Figure S1E). Moreover, Clock expression was significantly 
reduced (by about 25%) at CTO and CT4 in the small 
intestines of RORy~ ! ~ mice compared with WT controls 
(Supplementary Figure S1G). These results indicate that 
the regulation of Clock expression by RORs is tissue 
dependent and less stringent than that of Npas2. 
Examination of Clock expression in DKO mice showed 
that in contrast to single KO mice Clock expression was 
significantly reduced in the liver and BAT of DKO mice 
(Supplementary Figure SI I). 

Activation and in vivo ROR binding to the 
Npas2 promoter 

RORs bind as a monomer to ROREs consisting of the 
consensus sequence AGGTCA preceded by a 6-bp A/T 
rich region, in the promoter regulatory region of target 



genes (1-3). The murine Npas2 proximal promoter 
region contains 2 successive putative ROREs at —1372 
and —1399 that are spanned by 21 bp (Figure 2A). These 
ROREs are conserved between mouse and human Npas2. 
Although EMSA showed that RORs can bind these 
ROREs (27), the functionality of these ROREs in the 
intact promoter and whether RORs are recruited to this 
Npas2 promoter region in vivo have not been studied. To 
examine whether RORs were able to activate the Npas2 
proximal promoter, reporter gene assays were performed 
in human hepatoma Huh-7 cells and human kidney epi- 
thelial HEK293T cells with a pGL4.10-A^a?2(-1534/+81) 
reporter plasmid in which the Luc reporter was placed 
under the control of the — 1534/+81 proximal promoter 
region of Npas2. Figure 2B shows that in Huh-7 cells 
RORy was able to activate the Npas2 promoter more 
than 20-fold compared with transfection with the 
p3xFlag-CMV10 empty vector. Expression of RORa 
also greatly increased the activation of the Npas2{-\ 534/ 
+81) promoter. To determine which of the two sites was 
important in this activation, the effect of AGGTCA to A 
AATCA mutations in the RORE1 and RORE2 on the 
activation of the Npas2 promoter was examined. 
Mutation of RORE1 modestly reduced the activation of 
the Npas2 promoter by either RORa or RORy and 
mutation of the RORE2 decreased the activation more 
than the RORE1 mutation. Mutation of both ROREs 
almost completely abolished the activation of the 
A / /?av2(-1534/+81) promoter by RORs. Similar results as 
for Huh-7 were obtained in HEK293T cells (Figure 2C). 
During submission of this manuscript a different study 
provided evidence for a role of another RORE-like site 
(RORE3) at —1207 of the mouse Npas2 promoter in its 
activation by RORa (53). However, mutation of this 
RORE within the 7V>av2(-1534/+81) promoter had little 
effect on its activation by RORa or RORy and did not 
affect its repression by Rev-Erba, suggesting that this site 
plays a minor role in the regulation of Npas2 by these 
receptors (Supplementary Figure S2). Our observations 
indicate that RORE1 and RORE2, rather than RORE3, 
contributed to the activation of the Npas2(-1534/+&1) 
promoter by RORa and RORy, with a more prominent 
role for RORE2 than RORE1. 

Inhibition of A^>a,s2(-1534/+81) promoter activation by 
Rev-Erba and T0901317 

Earlier studies have shown that the nuclear receptor, 
Rev-Erba, which also binds ROREs, can compete with 
RORs for RORE binding (22,23,37,54). Figure 3A shows 
that co-expression of Rev-Erba repressed the activation 
of the Npas2 promoter by either RORa or RORy in a 
dose-dependent manner. Moreover, expression of Rev- 
Erba in Huh-7 cells repressed the endogenous activation 
of the Npas2(-l534/+8l) promoter (Figure 3B). 
Mutation of RORE1/2 abolished this repression 
indicating that this repression was mediated through 
these ROREs. These observations are consistent with 
the hypothesis that Rev-Erba functions as a repressor 
of Npas2 expression and competes with RORs for 
RORE binding. Recently, T0901317 was identified as 
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Figure 2. Regulation of Npos2 expression by ROR is mediated through ROREs. (A) Schematic of the location of the RORE1 and RORE2 in the 
proximal Npas2 promoter region. The bold and underlined sequences of these ROREs are conserved between the human and mouse genomes. The 
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an inverse agonist of RORa and RORy (7). It was 
therefore of interest to examine the effect of T0901317 
on the expression of Npas2 induction. As shown in 
Figure 3C, T0901317 was able to repress the activation 
of the Npasl promoter by RORy and RORa in Huh-7 
cells. 

Recruitment of RORs to the Npas2 promoter 

To determine whether RORy was recruited to the ROREs 
in the /V>aj2(-1534/+81) promoter in vivo, ChlP-qPCR 
analysis was carried out using liver isolated from WT 
and ROR knockout mice. Tissues were collected at 
CT22, a time close to the peak expression of RORy and 
Npasl mRNAs, and at CT10, a time at which both RORy 
and Npasl mRNAs are expressed at their lowest level. On 
this basis, one would expect higher recruitment of RORy 
to the Npasl promoter at CT22 than at CT10. The data in 
Figure 4A confirmed that RORy was more efficiently re- 
cruited to the —1380 region of the Npasl promoter in WT 
liver collected at CT22 compared with the liver collected at 
CT10. Primers targeting a region within the Gapdh gene 
were used as negative control and did not show any 
specific recruitment of RORy. As expected, at CT22, the 
level of RORy associated with the Npasl promoter in WT 
liver was significantly higher than in RORy~ ! ~ liver, while 
RORy was not recruited to a distal, non RORE- 
containing region of the Npasl gene. Although the loss 
of RORa had only a small effect on Npasl expression, 



ChIP analysis indicated that RORa was associated with 
the Npasl promoter. 

Npas2 activation and repression by RORs and Rev-Erba, 
respectively, in Hepal-6 cells 

To obtain further evidence for the involvement of RORa/y 
and Rev-Erba in the regulation of Npasl, we established 
murine hepatoma Hepal-6 cell lines stably expressing Flag- 
tagged RORa, RORy or Rev-Erba. Immunocytochemistry 
using anti-Flag M2 antibody showed that these receptors 
were expressed in these cells and localized to the nucleus 
(data not shown). Figure 5A shows that Npasl expression 
was dramatically induced in Hepal-6 cells expressing either 
RORa or RORy. In contrast, Rev-Erba overexpression 
significantly downregulated the expression of Npasl. 

To analyze the recruitment of RORa, RORy and 
Rev-Erba to the ROREs within the /V>as2(-1534/+81) 
promoter, we performed ChIP assay using the Hepal-6 
stable cell lines and an anti-Flag M2 antibody 
(Figure 5B). Our data showed that all three receptors 
were recruited to the —1380 region of the Npasl 
promoter, consistent with the conclusion that Npasl is 
directly regulated by RORs and Rev-Erba. 

Activation of the Npas2 promoter in relation to changes in 
DNA accessibility and histone H3K9 acetylation 

The FAIRE assay has been used to detect active regula- 
tory regions (50) and although the roles of histone 
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Figure 3. Rev-Erbtx and the ROR-inverse agonist T0901317 inhibited ROR-induced activation of the A^as2(-1534/+81) promoter. (A) Rev-Erba 
expression represses the activation of the Npas2 promoter by RORa and RORy, Huh-7 cells were transfected with p3xFlag-CMV10-RORy or 
p3xFlag-CMV10-RORa, pGL4.10-7V/«z.v2(-1534/+81) and increasing concentrations of p3xFlag-CMV10-Rev-Erba and 24h later were assayed for 
reporter activity. (B) The downregulation of basal Npas2 promoter activity by Rev-Erba was abrogated by mutations in the ROREs. Huh-7 cells 
were transfected with p3xFlag-CMV10-Rev-Erba and pGL4.10 driven by the WT 7V/?o.?2(-1534/+81) promoter or the promoter with the indicated 
RORE mutations. About 24 h later cells were assayed for reporter activity. (C) The inverse agonist, T0901317, represses the activation of the Npas2 
promoter by both RORy and RORa in Huh-7 cells. Data present mean ± SEM. 
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Figure 4. RORy and RORa are recruited to the Npas2 promoter. (A) Orcadian time-dependent recruitment of RORy to the Npas2 promoter. ChIP 
analysis was performed using an anti-RORy antibody and liver tissues (n = 4) isolated from WT and RORy~l~ mice at CT10 (low expression of 
RORy) and CT22 (high expression of RORy). Amplification of the Gapdh gene was used as a negative control. Data present mean ± SEM, 
***/>< 0.001 by ANOVA. (B) ChIP analysis was performed with anti-ROR antibodies using liver tissues (n = 4) isolated from WT, RORa ssjsg 
and RORy~l~ mice at CT22. QPCR amplification of a non RORE-containing distal site of the Npas2 gene was used as a negative control. Data 
present mean ± SEM, **P<0.01, ***/>< 0.001. 



acetylation, including H3K9Ac, are far from being 
understood levels of acetylation correlate positively 
with actively transcribed genes (55). To study the role 
of chromatin modulation in Npas2 regulation, we first 
examined the FAIRE and the H3K9Ac signals within 
the proximal Npas2 promoter region in DNA samples 
from liver isolated at CT10 and CT22. As shown in 
Figure 6, the FAIRE and H3K9Ac signals were 
markedly increased in DNA samples from liver isolated 



at CT22, a time at which Npas2 is most highly expressed, 
compared with that of CT10. These results are consistent 
with the concept that at CT22 this region of the Npas2 
promoter is accessible and associated with functionally 
active enhancers. Loss of either RORa or RORy expres- 
sion led to a reduction in FAIRE signal and H3K9 
acetylation at CT22 suggesting that ROR binding 
promotes open chromatin structure at the site of the 
Npas2 promoter. 



4776 Nucleic Acids Research, 2011, Vol. 39, No. 11 



A 180 




FAIRE 
Npas2 RORE1,2 



FAIRE 

Npas2 distal 




■ : A/pas2 RORE 
| : Npas2 distal 



Figure 5. Overexpression of RORot/y or Rev-Erba in Hepal-6 cells, 
respectively, induced or repressed Npas2 activation. (A) Npas2 gene 
expression was examined by QPCR analysis in Hepal-6 cells (n = 5) 
stably expressing empty vector, Flag-RORa, Flag-RORy or Flag-Rev- 
Erboc. The expression of Npas2 in Hepal-6(Empty) was normalized to 
1. Data present mean ± S EM, ***/>< 0.001 by ANOVA. (B) RORs 
and Rev-Erbc* were recruited to the Npas2 promoter in Hepal-6 
cells. ChIP analysis was performed with the Hepal-6 stable cell lines 
and anti-Flag M2 antibody. Hepal-6(Empty) served as a negative 
control. Data present mean ± SEM, ***/"< 0.001. 

Regulation of RORy and Npas2 expression in brown 
adipocytes 

Because of the observations that the expression of Npas2 
was most significantly affected in BAT of RORy~ ! ~ mice 
(Figure 1C) and that RORy is highly expressed in BAT 
compared with WAT and liver (Figure 7A), we were inter- 
ested in exploring the relationship between the expression 
of RORs and Npas2, and BAT differentiation and 
function. Analysis of RORy and Npas2 expression 
showed that both were induced about 80- and 50-fold, 
respectively, during differentiation of primary brown 
preadipocytes (Figure 7B and C). In contrast to RORy, 
Npas2 expression was repressed at later stages of differen- 
tiation. The induction of Npas2 mRNA expression was 
greatly reduced during differentiation of brown 
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Figure 6. Activation of the Npas2 promoter in relation to changes in 
DNA accessibility and histone H3K9 acetylation. (A) Accessibility to 
Npas2(ROREl,2), as assessed by FAIRE-QPCR analysis, correlates 
with Npas2 promoter activation. FAIRE-QPCR at 7V/m.s2(RORE1,2) 
and at distal downstream site of the Npas2 gene was performed using 
DNA samples isolated from liver of WT, RORa sg/sg and RORy- 1 - mice 
at CT10 and CT22. (B) Activation of the Npas2 promoter correlates 
with H3K9 acetylation. ChlP-QPCR was performed at the proximal 
Npas2 promoter (Npas2 -896/-813) and the distal site using 
anti-H3K9Ac antibody and ChlPed-DNA samples isolated from liver 
of WT, RORa sglsg and RORy- 1 - mice at CT10 and CT22. Data present 
mean ± SEM, *.P<0.05, ***P< 0.001. 

preadipocytes isolated from RORy~ ! ~ mice. However, 
the expression of the adipocyte marker genes, Glut4 and 
Fabp4, was at most time points not significantly different 
between WT and RORy~'~ adipocytes (Figure 7D and E) 
indicating that the downregulation of Npas2 in RORy~'~ 
brown adipocytes was rather specific. BAT plays an im- 
portant role in heat generation and a number of genes, 
including Mfsd2a and Dio2, are greatly induced during 
cold-induced thermogenesis (49,56). Figure 7F shows 
that the level of Npas2 mRNA in the BAT was increased 
about 7-fold during cold-induced thermogenesis and that 
this induction was inhibited by 90% in BAT of RORy~'~ 
mice. Loss of RORy also repressed the induction of Dio2 
and Mfsd2a during cold-induced thermogenesis. These 
data suggest that, in addition to its circadian regulation 
of Npas2, RORy plays a critical role in the regulation of 
Npas2 expression during BAT differentiation and 
cold-induced thermogenesis. 

To study the role of RORy in the regulation of Npas2 
in brown adipocytes further, we examined its expression 
in immortalized brown preadipocyte cell lines stably 
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Figure 7. The induction of Npas2 expression during brown adipocyte differentiation was regulated by RORy. (A) Comparison of RORa and RORy 
expression in the WAT, BAT and several other tissues. RORa and RORy were highly expressed in the WAT and BAT, respectively. (B-E) RORy 
mRNA was highly induced during brown adipocyte differentiation. Primary cultures of BAT preadipocytes obtained from WT and RORy~^~ mice 
were induced to differentiate as described in 'Materials and Methods' section. RNA was isolated at the times indicated and the expression of RORy, 
Npas2 and BAT differentiation markers, such as Glut4 and Fabp4, was examined by QPCR analysis. The data are an average of three isolates of 
brown adipocytes from each WT and RORy~ f ~ BAT, respectively, e.g. preadipocytes in exponential phase. (F) The induction of Npas2 expression 
during cold-induced thermogenesis was inhibited in the BAT of RORy~ l ~mice. WT and RORy~l~ mice (n = 5) were placed at 4°C for 4h before 
BAT was collected for analysis at CT18. Data present mean ± SEM, **P<0.01, ***p< 0.001 by ANOVA. 



expressing empty vector, RORy or the mutant 
RORy(E502Q), referred to as BAT(E), BAT(RORy) and 
BAT(RORyE502Q). As shown in Figure 8A, Npas2 was 
induced at 6- to 7-fold higher levels during differentiation 
of BAT(RORy) cells compared with BAT(E) cells. In 
contrast, Npas2 was not elevated in BAT(RORyE502Q) 
cells expressing an RORy mutant lacking transactivation 
function. Treatment of BAT(RORy) cells with the inverse 
agonist, T0901317, significantly reduced the induction of 
Npas2 mRNA (Figure 8A). 

ChIP analysis with anti-Flag M2 antibody and 
BAT(RORy) or BAT(E) cells showed strong recruitment 
of Flag-RORy to the RORE-containing region of the 
Npas2{-\51>AI+%\) promoter (Figure 8B). These data 
support our conclusion that RORy is recruited to the 
Npas2{-\51>AI+%\) promoter. The inverse agonist, 
T0901317, did not significantly influence the recruitment 
of RORy to the Npas2 promoter suggesting that 
T0901317-bound RORy was still recruited to ROREs 
within the Npas2 promoter. 



DISCUSSION 

In this study, we provide evidence that RORs, RORy in 
particular, are important regulators of Npas2 expression. 



In WT mice, Npas2 displayed a strong oscillatory expres- 
sion pattern in several tissues, including the kidney, liver 
and BAT, with peak expression at CT20-CT24/0 
(Figure 1) (24,51). Although the rhythmicity of Npas2 
expression in these tissues was maintained in RORy~ l ~ 
mice, the peak level of Npas2 mRNA expression was sig- 
nificantly reduced. Except for the kidney, loss of RORa 
had relatively little effect on Npas2 expression in the 
tissues analyzed despite its recruitment to the Npas2 
promoter. That loss of RORa has little effect on Npas2 
expression might be related to weak RORa activity and 
involve low levels of RORa protein bound relatively to 
RORy or posttranslational modifications of RORa. 
Additional evidence that RORs regulate Npas2 expression 
was provided by data showing that overexpression of 
either RORa or RORy in murine hepatoma Hepal-6 
cells greatly induced the expression of Npas2 mRNA 
(Figure 5A). 

We previously reported a certain degree of functional 
redundancy between the regulation of gene expression by 
RORa and RORy (16). For example, the hepatic expres- 
sion of several Phases I and II metabolic genes was con- 
siderably more affected in DKO than in single knockout 
mice. Analysis of Npas2 expression in DKO mice 
indicated no significant functional redundancy between 
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Figure 8. RORy-activated Npas2 gene expression in brown adipocytes 
through a direct mechanism. (A) Npas2 expression was examined in 
brown adipocyte cell lines, BAT(E), BAT(RORy) and 
BAT(RORyE502Q) stably expressing empty vector, RORy or the 
RORyE502Q mutant containing an inactive activation domain, respect- 
ively. Expression of Npas2 was compared between brown preadipocytes 
at the exponential phase (U) and at differentiation Day 5 (D) as 
indicated. BAT(RORy) cells were also treated at Day 4 for 24h with 
10 uM of the inverse agonist T0901317. (B) ChIP analysis was per- 
formed using anti-Flag M2 antibody, BAT(E) cells (E) and 
BAT(RORy) cells (R) treated with or without T0901317 as described 
under (A). BAT(E) cells and the distal region were used as a negative 
controls. The experiment was performed in triplicate. Data present 
mean ± SEM, ***/>< 0.001 by ANOVA. 

the RORa and RORy in regulating Npas2. These data are 
consistent with our conclusion that in these tissues Npas2 
expression is regulated selectively by RORy with a minor 
role for RORa. In contrast to Npas2, the expression of 
Clock was significantly reduced in the liver and BAT of 



DKO mice (Supplementary Figure SI I), while the expres- 
sion was not changed in single knockout mice suggesting a 
redundant function of RORa and RORy in regulating 
Clock. 

RORs can regulate gene transcription directly by 
binding ROREs in the regulatory region of target genes 
or through indirect mechanisms (1). The 1534-nt proximal 
promoter region of the Npas2 gene contains three putative 
ROREs (27). In this study, we demonstrate that Npas2 is a 
direct ROR target gene and that its transcriptional acti- 
vation involves the recruitment of RORs to RORE1 and 
RORE2. This is supported by reporter gene analysis 
showing that expression of either RORa and RORy 
activated the Npas2(- 1534/ '+81) proximal promoter. 
Mutation of either RORE1 or RORE2 partially reduced 
the Npas2 promoter activity, whereas mutation of both 
ROREs almost totally abolished activation of the Luc 
reporter suggesting that optimal Npas2 induction by 
RORy requires both ROREs. In contrast to an earlier 
report, mutations in RORE3 had little effect on ROR- 
induced activation of the Ayja.v2(-1534/+81) promoter 
suggesting that it plays a minor role in the regulation of 
Npas2 by RORs. We further demonstrated by ChIP 
analysis that RORa and RORy were recruited to these 
ROREs in the Npas2 promoter in WT liver, Hepal-6 
and BAT cells. As expected no significant recruitment of 
RORa or RORy was observed in RORa s/ss or RORy~'~ 
liver, respectively. Our data further showed that the re- 
cruitment of RORy was CT dependent. In contrast to 
CT22, when RORy is expressed at a high level, RORy 
was not found to be associated with the Npas2 promoter 
at CT10, when RORy expression is low. Collectively, these 
data suggest that the regulation of Npas2 transcription by 
RORy occurs by a direct mechanism. Transcriptional 
activation is associated with increased chromatin accessi- 
bility, as determined by FAIRE, and histone acetylation, 
including H3K9 acetylation (50,55). Our data showed that 
the FAIRE signal and H3K9 acetylation are significantly 
increased at the Npas2 promoter in liver at CT22 
compared with CT10 consistent with increased accessibil- 
ity and the high level of transcriptional activity of the 
Npas2 promoter at CT22. Loss of RORy significantly 
reduced the FAIRE signal and H3K9 acetylation at the 
Npas2 proximal promoter consistent with a relationship 
between RORy binding, open chromatin structure and 
transcriptional activation of Npas2. Despite the fact that 
loss of RORa has little effect on Npas2 expression in liver, 
it also caused a reduction in the FAIRE signal and H3K9 
acetylation suggesting a role for RORa plays a role in 
chromatin remodeling in Npas2. This is consistent with 
the ChIP analysis showing binding of RORa to the 
Npas2 ROREs (Figure 4) (53). Further studies are 
required to understand the role of RORa in Npas2 gene 
regulation. Chromatin remodeling might occur prior to 
and/or following ROR binding. The observation that 
loss of RORs did not totally inhibit Npas2 expression 
and reduced FAIRE signal and H3K9 acetylation by 
about 50% suggests that other transcription factors are 
involved in the positive regulation of Npas2 and chroma- 
tin remodeling of the Npas2 promoter. Thus, at CT22 
RORs may bind to this region when it is already quite 
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accessible through the activity of other transcription 
factors. RORs through their recruitment of histone 
acetyltransferases, such as CBP/p300, likely contribute 
to the increased histone acetylation and chromatin 
remodeling. 

Recent studies have demonstrated that RORs function 
as ligand-dependent transcription factors (5-8,58). 
T0901317 functions as an inverse agonist for both 
RORa and RORy and has been reported to effectively 
inhibit RORa- and RORy-mediated induction of Thl7 
differentiation and IL-17 expression. In this study, we 
demonstrate that T0901317 inhibited the activation of 
the Npas2 promoter by RORa and RORy in hepatoma 
Huh-7 cells as well as in brown adipocytes. ChlP analysis 
showed that T0901317 did not significantly affect the re- 
cruitment of RORy to the ROREs in the Npas2 proximal 
promoter. The recruitment of RORy to the Npas2 
promoter in the presence of T0901317 indicates that 
T0901317 does not affect the interaction of RORs with 
these ROREs. The inhibition of transcription is likely 
due to a T0901317-induced change in the conformation 
of the ROR receptors that affects their interaction with 
co-activators and subsequently their transcriptional 
activity. This hypothesis is supported by a recent study 
showing that T0901317 inhibits the interaction of RORs 
with the co-activator SRC-1 (7). 

The nuclear receptors, Rev-Erba and Rev-Erb(3, have 
been reported to act as transcriptional repressors (59). 
Moreover, it is well-established that Rev-Erb receptors 
bind response elements similar to RORE and can 
compete with RORs for RORE binding (22,23,37,54). In 
this study, we demonstrated that Rev-Erba can also 
function as a repressor of Npas2 expression. In contrast 
to the effects of RORa or RORy, overexpression of 
Rev-Erba in Hepal-6 cells caused a reduction in Npas2 
mRNA level (Figure 5A). In addition, overexpression of 
Rev-Erba inhibited the basal Npas2(—1534/+81) promoter 
activity (Figure 3B). Mutation of either RORE1 or 
RORE2, but not RORE3, reduced this inhibition, while 
the double mutation totally abrogated the repression by 
Rev-Erba, suggesting that the inhibitory effect of 
Rev-Erba is mediated through these two ROREs. This 
was supported by ChlP analysis showing recruitment of 
Rev-Erba to the RORE-containing region of the Npas2 
promoter (Figure 5B). Thus, as RORs, Rev-Erbs play a 
critical role in the regulation of Npas2 expression through 
recognition of the same ROREs in the Npas2 promoter. It 
is likely that competition between RORs and Rev-Erba 
for RORE1 and RORE2 is part of the mechanism that 
regulates the oscillatory expression of Npas2. This is 
corroborated by data showing that Rev-Erba effectively 
antagonized the activation of the Npas2 promoter by 
RORs (Figure 3A). 

We were interested in relating the regulation of Npas2 
by RORy and Rev-Erba to their rhythmic pattern of ex- 
pression. The oscillatory expression pattern of Npas2 and 
RORy are in almost opposite phase compared with that of 
Rev-Erba, while RORy reaches a peak slightly earlier than 
Npas2 (Figures 1 and 9) (16,20,24,51). The lowest level of 
Rev-Erba mRNA expression occurred at CT16-CT24 
when Npas2 and RORy are expressed at peak levels, 



whereas at CT8-CT10 the peak level of Rev-Erba expres- 
sion coincides with the lowest level of RORy and Npas2 
expression. These patterns of expression are consistent 
with the hypothesis that at CT8-10 downregulation of 
RORy and induction of Rev-Erba may act cooperatively 
in repressing Npas2 expression, while induction of RORy 
and repression of Rev-Erba jointly promote Npas2 expres- 
sion at CT16-24 (Figure 9). A recent study demonstrated 
that loss of Bmall abrogates Rev-Erba expression and as a 
result alleviates the repression of RORy expression by 
Rev-Erba (24), while RORa has been demonstrated to 
positively regulate Rev-Erba (22,60,61). Thus, RORy 
and Rev-Erba are part of a regulatory loop that regulates 
the circadian pattern of expression of Npas2 (Figure 9). 
The opposite phases of the circadian expression of these 
receptors and their competition for RORE binding col- 
lectively participate in the regulation of the oscillatory 
pattern of Npas2 expression. 

As its paralog Clock, Npas2 forms a heterodimeric 
complex with Bmall and regulates gene transcription by 
interacting with E-boxes in the regulatory region of target 
genes (39). Clock and Npas2 have been reported to exhibit 
overlapping roles in the regulation of the circadian clock 
genes in the SCN and circadian-controlled genes in liver 
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Figure 9. Schematic presentation of the molecular mechanism for the 
antagonistic regulation of Npas2 circadian expression by RORy and 
Rev-Erba. (Lower panel) Oscillatory pattern of expression of Npas2, 
RORy and Rev-Erba in liver. (Upper panel) RORy functions as a 
positive regulator of Npas2 transcription in vivo and in cultured cells 
by binding to two ROREs in its proximal promoter region. RORy is 
efficiently recruited to these ROREs at CT18-CT24/0 when RORy and 
Npas2 are highly expressed. Rev-Erba, which peak expression occurs 
around CT8 (37), represses Npas2 expression through competitive 
binding to the same ROREs. Rev-Erba has also been reported to 
repress RORy expression, while RORy can activate Rev-Erba expres- 
sion. The reduced RORy expression and increased Rev-Erba expression 
at CT8-10 might act synergistically to reduce Npas2 circadian expres- 
sion, while the inverse occurs at CT16-24. Thus, RORy and Rev-Erba 
are parts of a regulatory loop regulating each other and consequently 
the circadian expression of Npas2. Because loss of RORs does not 
totally block Npas2 expression, other transcription factors likely 
regulate its transcription as well. 
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(39-41). In addition to its role in circadian regulation, 
epidemiologic studies have provided evidence for a role 
of Npas2 in several other biological processes. SNPs in 
Npas2 have been linked to increased risk of cancer, 
unipolar major depression, metabolic syndrome and 
hypertension (42-47). Npas2 has been implicated in the 
regulation of the DNA damage response pathway and 
several cell cycle and DNA repair genes suggesting a 
role for Npas2 as a tumor suppressor (45,46). In this 
study, we demonstrate that both RORy and Npas2 are 
induced during differentiation of preadipocytes into 
brown fat adipocytes. In addition, we show that Npas2 
expression is enhanced in BAT during cold-induced 
thermogenesis, suggesting a regulatory role for Npas2 in 
this physiological process. However, Npas2 expression 
was strongly downregulated in RORy~ ! ~ brown adipo- 
cytes during differentiation and in BAT of RORy~ l ~ 
mice during cold-induced thermogenesis. In addition, we 
showed that T0901317 was able to inhibit the induction of 
Npas2 by RORy in the brown adipocytes. Collectively, 
these data support the conclusion that RORy functions 
as a transcriptional activator of Npas2 in BAT. 
Although our study suggests a role for Npas2 in thermo- 
genesis, future studies are required to determine the 
precise function of Npas2 and RORy in BAT cells. 
Taken the role of Rev-Erbs, RORs and BAT in metabolic 
syndrome, these data are of interest in the light of reports 
showing that Npas2 gene variants have been linked to 
increased risk of metabolic syndrome (16,17,22,43,62). 

Our study shows that RORy regulates Npas2 expression 
in vivo. Loss of RORy had little effect on the rhythmicity, 
but reduced the peak level of Npas2 expression. RORy 
regulates the transcription of Npas2 by binding to two 
ROREs in its proximal promoter region and enhances 
the accessibility to the proximal Npas2 promoter 
enhancer region. The activation of Npas2 promoter by 
RORs is repressed by co-expression with Rev-Erba or 
addition of the inverse agonist T0901317. These observa- 
tions suggest that RORs and Rev-Erba may be implicated 
in the modulation of Npas2-controlled physiological 
processes and Npas2-linked pathologies, including circa- 
dian behavior disorders, metabolic syndrome and tumori- 
genesis. Thus, changes in the level of the physiological 
ligands of RORs and Rev-Erba may play an important 
role in the development of these pathologies. In addition, 
because RORs function as ligand-dependent transcription 
factors, increased Npas2 expression by ROR agonists may 
provide a novel therapeutic approach for disease in which 
Npas2 has been implicated, including metabolic 
syndrome, hypertension and depression. 
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